clear,clc
x=100:100:500;
y=100:100:400;
z=[636    697    624    478   450  
   698    712    630    478   420
   680    674    598    412   400
   662    626    552    334   310];
pp=csape({x,y},z')
xi=100:10:500;yi=100:10:400;
cz=fnval(pp,{xi,yi});
[i,j]=find(cz==max(max(cz)))  %找最高点的地址
x=xi(i),y=yi(j),zmax=cz(i,j)  %求最高点的坐标 
